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ABSTRACT 

We present a quantitative estimate of the relative AGN/starburst content in a sample 
of 59 nearby [z < 0.15) infrared bright ULIRGs taken from the 1 Jy sample, based on 
infrared L-band (3-4 /im) spectra. By using diagnostic diagrams and a simple decon- 
volution model, we show that at least 60% of local ULIRGs contain an active nucleus, 
but the AGN contribution to the bolometric luminosity is relevant only in ~ 15 — 20% 
of the sources. Overall, ULIRGs appear to be powered by the starburst process, re- 
sponsible for > 85% of the observed infrared luminosity. The subsample of sources 
optically classified as LINERs (31 objects) shows a similar AGN/starburst distribu- 
tion as the whole sample, indicating a composite nature for this class of objects. We 
also show that a few ULIRGs, optically classified as starbursts, have L-band spectral 
features suggesting the presence of a buried AGN. 
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1 INTRODUCTION 

Ultraluminous Infrared Galaxies (ULIRGs, Lir ~ Lbol > 
10 12 Lq) have been studied at all wavelengths, from radio to 
hard X-rays, with the primary aim of estimating the relative 
contribution of accretion and star formation activity to the 
bolometric luminosity. 

Recently, L-band (3-4 pm) spectroscopy provided new 
powerful tools to perform this study. Imanishi & Dudley 
(2000), Risaliti et al. (2006, hereafter R06), Imanishi et 
al. (2006, hearafter 106), Imanishi et al. (2008), Sani et 
al. (2008) analyzed nearby bright ULIRGs from the IRAS 
Bright Galaxy Sample (Sanders et al. 2003) and the IRAS 
1 Jy sample (Kim & Sanders 1998) and demonstrated that 
the AGN and starburst components can be disentangled in 
L-band spectra, using several indicators: 

- 3.3 jUm emission feature: this emission is due to Pol- 
icyclic Aromatic Hydrocarbon (PAH) molecules, and is 
prominent in starburst-dominated sources (equivalent width 
EW~ 110 nm), while it is weak or absent in AGN-dominated 
sources. This observational result is motivated by the differ- 
ent intrinsic emission of AGNs and starbursts: the intense 
X-ray radiation of the former is believed to destroy the PAH 
molecules, which can instead survive in the starburst radia- 
tion field. 

- 3.4 fim absorption feature: in ULIRGs, the presence of 
this feature with an optical depth T3.4 > 0.2 is observa- 
tionally associated to heavily obscured AGNs. This is due 



to the centrally concentrated emission of AGNs, which can 
be effectively covered by large columns of dust. Imanishi & 
Maloney (2003) demonstrated that a large-scale source, such 
as a starburst, interspersed within the dust responsible for 
absorption, cannot have a 3.4 /im absorption feature with 
r > 0.2. 

- Continuum slope: ULIRGs containing heavily obscured 
AGNs show a much redder L-band continuum than unob- 
scured AGNs and starbursts. This is due to the reddening 
effect of the covering dust. If the continuum is modeled with 
a simple power law in the A- fx plane, i.e. f\(x\ F , the typi- 
cal values for the slope P are F^-0.5 for unobscured AGNs 
and pure starburst, while F>1 for heavily obscured AGNs 
(R06). 

- Bolometric ratios: the intrinsic ratio between the 3 /im and 
bolometric luminosity is ~ 100 times higher in pure AGNs 
than in starbursts (R06). The measured ratio in ULIRGs is 
therefore in itself an indicator of the possible presence of an 
AGN component. 

The power of L-band spectroscopic analysis lays in the capa- 
bility of detecting heavily obscured AGNs, which are missed 
at other wavelenghs. The combination of the properties 
listed above implies that an AGN with an L-band optical 
depth Al ~ 1 — 5 (corresponding to 71^^20-100, i.e. com- 
pletely absorbed from the UV to the near-IR) can still be 
detectable in the L-band spectrum of ULIRGs, even when 
its contribution to the bolometric luminosity is negligible. 
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A quantitative application of the above diagnostics has been 
proposed by R06, who assumed that the L-band spectra can 
be reproduced as the combination of a fixed starburst tem- 
plate and an AGN template with variable dust absorption. 
As a consequence, the model has two free parameters: the 
fraction of the AGN contribution to the L-band luminosity, 
a, and the optical depth tl of the dust absorbing the AGN 
component. These two parameters can be estimated from 
the observed continuum slope, F, and the equivalent width 
EW3.3 of the 3.3 fim PAH emission feature. 

In this paper we present the application of this method 
to a complete sample of 59 ULIRG with available L-band 
spectroscopic observations. The sample consists of all the 
ULIRGs in the IRAS 1 Jy sample with redshift 2 < 0.15. 
The optical classification (Veilleux et al. 1999) is AGN for 
14 sources (25%), starburst for 12 sources (22%) and LINER 
for 31 sources (56%). Two sources are not optically classi- 
fied. Our work is mostly based on published data, which 
already provided an AGN/starburst classification, based on 
the diagnostic features summarized above. The significant 
improvement presented here consists in moving from a sim- 
ple ULIRGs classification as starburst or AGN, to a quan- 
titative estimate of the contribution of the two components 
to the bolometric luminosity. 

A similar approach has also been successfully used by 
Nardini et al. (2008) on Spitzer/IRS 5-8 fim spectra of a 
sample of local ULIRGs largely overlapping with the one 
presented here. We will discuss our results in comparison 
with those of Nardini et al. 2008 in Section 3. 



2 DATA ANALYSIS 

The observations in the northern hemisphere have been per- 
formed with the SUBARU (45 sources) and IRTF (6 sources) 
telescopes, and have been presented in 106. In the southern 
hemisphere the observations were made with the VLT (eight 
sources). Four of these sources have been presented by R06 
and Risaliti et al. 2006B, while the remaining four were ob- 
served in early 2006 with ISAAC, with the same configura- 
tion as for the first eight sources, and reduced following the 
steps described in R06. We refer to R06 and 106 for a de- 
tailed description of the data reduction and a visual analysis 
of the spectra. 

We performed a homogeneous analysis of all the sources 
in our sample, fitting each fx spectrum with a standard 
X 2 minimization using stated observational uncertainties. 
The adopted model consists of a power law continuum, and 
a Gaussian emission feature at rest frame wavelength of 
3.3 ^m. For each source, the spectral interval used for the 
fit is A = 3.25 x (1 + z) — 4.1 /im, i.e. the region containing 
the 3.3 fim emission feature and the continuum at longer 
wavelengths up to the end of the L-band. We excluded the 
wavelengths shortward of the PAH emission feature to avoid 
regions of bad atmospheric transmission, and possible dis- 
tortions of the continuum slope due to ice absorption fea- 
tures (see R06 and 106 for more details). In a few spectra an 
hydrocarbon absorption feature at rest frame wavelength 
A res t = 3.4 — 3.5 fim is clearly present. In these cases we 
added a Gaussian absorption component in our fits. In all 
the other cases we checked that this extra component is 
not required statistically. In Fig. 1 we show a few examples 
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Figure 1. Six examples of our data and best fit models, including 
objects among the highest (B) and lowest (F) S/N. A, B and C: 
sources intrinsically dominated by a heavily obscured AGN. The 
starburst component is however strong in the observed spectra in 
A and C. Note that C is optically classified as starburst. D and E: 
objects containing a weakely obscured AGN, which is dominant 
in the L-band, but not in the Bolometric luminosity (See Table 1). 
F: Starburst-dominated object. 



of the data and spectral fitting for representative objects 
in our sample. In Table 1 we list the best fit parameters 
r and EW3.3 for each source. A third important parame- 
ter which can be directly inferred from the analysis is the 
ratio between the luminosity in the L-band and the bolo- 
metric one, obtained from IRAS fluxes. As shown in R06, 
these three parameters are enough to provide a classifica- 
tion of the sources with respect to their AGN/starburst con- 
tent. This is illustrated in Fig. 2 and 3. In Fig. 2 we show 
the position of our sources in the T-EW3.3 plane. We also 
show the regions occupied by obscured AGN, unobscured 
AGN, and starbursts, based on the calibration done in R06. 
Fig. 3 shows the measured L-band to bolometric flux ratio, 
■Rl=A/a(3.3 ^m)/F(8-1000 (jxa) versus EWz.z. We note that 
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Source (instrument a ) 


r 


EW 


T 




a BO L 


T 3.4 


Cl(L) b 


Cl(Sp) c 


IRAS 00188-0856 (SJ 


6.2+0.1 


76+ 13 


7.7+ 1.4 


0.995+ 0.003 


0.91+ 0.06 


0.36+0.06 


AGN* 


AGN 


IRAS 02411+0353E (V) 


-2.4+2.2 


177+ 51 


0° 


< 0.1" 


< 10~ 3e 


A f 

J 


~ 


SB 


IRAS 02411+0353W (VJ 


-2.5+3.0 


183+ 80 


< 3 


< 0.96 


< 0.18 


A f 

J 




SB 


IRAS 03250+1606 (S) 


-1.8+1.1 


83+ 22 


Q d 


<0.52 


<0.01 


A f 


AGN* 


SB 


IRAS 04103-2838 (V) 


-0.1+0.2 


59+ 6 


0.7+ 0.4 


0.63+0.07 


0.02+0.005 


A f 

J 




AGN* 


IRAS 08572+3915 (I) 


1.9+0.1 


1+ 1 


1.7+ 0.1 


> 0.99 


> 0.71 


0.88+0.05 


AGN 


AGN 


IRAS 09039+0503 (S) 


0.3+0.8 


257+ 88 


e 


< 0.1 e 


< 10 — 3e 


J 


AGN* 


AGN* 


IRAS 09116+0334 (S) 


-1.9+0.3 


77+ 18 


Q d 


0.3+0.2 


n AA+0.015 

a02 It>.015 


J 


AGN* 


SB 


IRAS 09539+0857 (S) 


-0.5+0.2 


101+ 21 


<0.2 


<0.09 


<0.01 


f 


SB 


AGN* 


IRAS 10378+1108 (S) 


-2.8+1.2 


43+ 24 


d 


0.6+0.2 


n ao+o.02 


/ 


SB 


AGN* 


IRAS 10485-1447 (S) 


3.4+0.2 


41+ 3 


3.6+ 0.3 


0.98+ 0.01 


0.37+ 0.03 


0.25+0.08 


AGN* 


AGN* 


IRAS 10494+4424 (S) 


-0.5+0.2 


126+ 15 


<0.2 


<0.11 


<0.01 


0^ 


AGN* 


SB 


IRAS 11095-0238 (S) 


1.5+0.1 


127+ 26 


5.7+ 3.9 


>0.95 


> 0.20 


of 


SB 


AGN 


IRAS 12112+0305NE (V) 


0.1 + 1.0 


69+ 30 


< 3 


0.63+ 0.35 


< 0.04 


o' 


SB 


SB 


IRAS 12112+0305SW (V) 


-0.9+0.7 


< 300 


- 


- 


- 


of 


- 


- 


IRAS 12359-0725 (S) 


2.3+0.3 


83+ 17 


4.6+ 1.8 


0.97+ 0.02 


0.25+ 0.15 


of 


AGN* 


AGN* 


IRAS 12127-1412 (S) 


2.1+0.1 


21+ 6 


2.0+ 0.1 


0.97+ 0.01 


0.25+ 0.06 


of 


AGN 


AGN 


IRAS 13335-2612 (V) 


2.0+1.5 


220+ 146 


e 


< 0.1" 


< 10~ 3e 


of 


- 


SB 


IRAS 14252-1550E (S) 


0.1 + 1.6 


85+ 32 


- 






of 
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IRAb 14252-looOW (b) 


0.1+0.3 


91+ 20 


<1.8 


A r r _l_ A OA 
U.OO+ 0.30 


A A1 J_ A AAO 
0.01+ 0.008 
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AGIN 


SB 


IRAS 14348-1447 (V) 


0.4+0.4 


110+ 30 


1.2+ 1.1 


<0.72 


<0.03 


of 


AGN* 


AGN* 


Arp 220 (I) 


0.4+0.2 


79+ 9 


1.8+ 0.7 


0.69+ 0.10 


0.02+ 0.01 


of 


SB 


AGN* 


IRAS 16090-0139 (S) 


2.5+5.8 


< 180 


_ 


_ 


_ 


of 


AGN* 


AGN* 


IRAS 16468+5200 (S) 


-2.8+1.6 


184+ 70 


_ 


< 0.54 


< 0.01 


of 


SB 


AGN* 


IRAS 16487+5447 (S) 


-2.0+1.4 


126+ 37 


_ 


_ 


_ 


of 


AGN* 


SB 


IRAS 17028+5817 (S) 


-2.1+0.2 


83+ 11 




0.24+0.10 


0.003+0.001 


of 


AGN* 


SB 


IRAS 17044+6720 (S) 


1.3+0.1 


18+ 4 


1.4+ 0.1 


0.95+ 0.01 


0.17+ 0.03 


0.14+0.09 


AGN 


AGN 


IRAS 21219-1757 (S) 


-0.7+0.1 


1+ 0.5 


d 


0.99+ 0.005 


0.47+ 0.10 


0' 


_ 


_ 


IRAS 21329-2346 (S) 


-0.4+0.7 


79+ 16 


< 1.5 


< 0.74 


< 0.02 


0^ 


AGN* 


AGN* 


IRAS 23234+0946 (S) 


0.6+0.1 


146+ 14 


0" 


< 0.1" 


< 10~ 3e 


0' 


SB 


SB 


IRAS 23327+2913 (S) 


0.1+0.8 


52+ 16 


< 2.0 


0.70+ 0.20 


0.023+0.017 


of 


SB 


AGN* 




-0.5+0.1 


82+ 3 


<0.2 


A OfiX A nti 
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A nAQ-Ln AA1 

u.uuo+u.uul 


n/ 


SB 


SB 


IRAS 10190+1322W(S) 


-2.2+0.4 


79+ 20 


d 


< 0.46 


< 0.02 


o' 


SB 


SB 


IRAS 11387+4116 (S) 


-1.1+0.4 


106+ 31 


Q d 


<0.45 


<0.02 


o/ 


SB 


SB 


IRAS 11506+1331 (S) 


2.5+0.2 


87+ 10 


5.4+ 1.3 


0.98+ 0.01 


0.35+ 0.12 


0^ 


AGN 


AGN 


IRAS 13509+0442 (S) 


0.2+1.3 


157+ 35 


< 3.4 


<0.97 


<0.28 


of 


SB 


SB 


IRAS 13539+2920 (S) 


-0.6+0.6 


94+ 20 


< 2.8 


<0.73 


<0.03 


of 


SB 


SB 


IRAS 14060+2919 (S) 


-0.2+0.2 


185+ 15 


C 


< 0.1" 


< 10~ 3e 


of 


SB 


SB 


IRAS 15206+3342 (S) 


0.2+0.2 


71+ 14 


1.1+ 0.8 


0.64+ 0.15 


0.02+ 0.01 


of 


SB 


SB 


IRAS 15225+2350 (S) 


1.3+0.6 


28+ 13 


1.5+ 0.8 


0.93+ 0.04 


0.12+ 0.07 


of 


AGN* 


AGN* 


IRAS 16474+3430 (S) 


0.5+0.4 


103+ 29 


4.8+ 13.4 


0.90+ 0.48 


0.08+ 0.36 


of 


AGN* 


SB 


IRAS 20414-1651 (S) 


-1.4+1.0 


77+ 18 


1.5+ 3.4 


<0.66 


<0.02 


of 


SB 


SB 


IRAS 21208-0519 (S) 


-0.4+3.1 


169+ 102 








of 


SB 


SB 


IRAS 05189-2524 (V) 


0.1+0.1 


7+ 1 


0.4+ 0.1 


0.96+ 0.01 


0.18+ 0.01 


of 


AGN 


AGN 


IRAS 08559+1053 (S) 


-0.2+0.1 


13+ 4 


0.3+ 0.1 


0.91+ 0.03 


0.09+ 0.03 


of 


AGN 


- 


IRAS 12072-0444 (S) 


1.2+0.2 


86+ 19 


3.5+ 1.9 


0.90+ 0.08 


0.08+ 0.06 


0.57±0.11 


AGN 




Mrk 273 (I) 


0.1+0.1 


31+ 5 


0.5+ 0.2 


0.81+ 0.04 


0.04+ 0.01 


o' 


AGN 


AGN 


IRAS 13443+0802 (S) 


-0.6+1.1 


111+ 56 


< 2.0 


< 0.5 


< 0.01 


o/ 


SB 


_ 


PKS 1345+12 (I) 


0.7+0.1 


7+ 3 


0.9+ 0.1 


0.97+ 0.01 


0.25+ 0.08 


o/ 


AGN 


AGN 


IRAS 15130-1958 (S) 


0.2+0.1 


8+ 1 


0.5+ 0.1 


0.95+ 0.01 


0.17+ 0.05 


o' 


AGN 




IRAS 17179+5444 (S) 


1.3+0.3 


24+ 13 


1.5+ 0.4 


0.94+ 0.04 


0.14+ 0.09 


1 + 0.5 


AGN 




Mrk 1014 (I) 


-0.7+0.1 


8+ 1 


d 


0.91+ 0.01 


0.09+ 0.01 


0^ 


AGN 




IRAS 07599+6508 (S) 


-0.3+0.1 


1+ 1 


0.2+ 0.05 


> 0.99 


> 0.7 


o' 


AGN 




IRAS 11598-0112 (S) 


-0.3+0.2 


5+ 3 


0.1+ 0.1 


0.96+ 0.03 


0.19+ 0.10 


o' 


AGN 




Mrk 231 (I) 


-0.4+0.1 


3+ 0.5 


0.1+ 0.05 


0.98+ 0.01 


0.28+ 0.02 


o' 


AGN 


AGN 


IRAS 15462-0450 (S) 


0.1+0.2 


15+ 5 


0.4+ 0.2 


0.90+ 0.04 


0.09+ 0.03 


o' 


AGN 




IRAS 12127-1412 (S) 


2.0+0.1 


21+ 6 


2.0+ 0.1 


0.97+ 0.01 


0.25+ 0.06 


0.25+0.07 


AGN 


AGN 


IRAS 14197+0813 (S) 


-0.6+0.5 


130+ 16 


<0.9 


<0.67 


<0.02 


o/ 






IRAS 14485-2434 (S) 


0.7+0.1 


65+ 8 


1.7+ 0.3 


0.80+ 0.04 


0.04+ 0.01 


of 


AGN 





Table 1. Results from AGN-SB decomposition as described in the text. The horizontal lines divide the sources in groups according to 
their optical classification: LINERs (first group), starbursts (second group), AGN (third group) and unclassified (the last two sources). 
The error on oibol does not take into account the uncertainties on the L-band to bolometric ratios. A more realistic estimate can be 
estimate from the scatter in Fig. 5B (see text for more details). a : Telescope/instrument of the observation: S: Subaru; V: VLT-ISAAC; I: 
IRTF. b : Classification from "standard" L-band diagnostics, (from 106), based on the PAH EW and the presence of absorption features. 
An asterisk indicates a weak evidence for the presence of the AGN. c : Classification based on Spitzer-YRS spectroscopy (from Imanishi ct 
al. 2007). d : Absorption is assumed to be r = 0, due to the steep continuum (r < —0.5 at >90% confidence level). e : Values and upper 
limits obtained assuming the source is pure starburst, due to the high equivalent width of the 3.3/im feature (EW>110 nm at > 90% 
confidence level), f; 3.4 /im absorption not required to adequately fit the spectrum. 
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Rl is in principle a rather powerful indicator, given the huge 
difference between the L-band to bolometric ratios for pure 
AGNs and pure starbursts (of a factor of about 100, R06). 
However it is not capable of giving an absolute estimate 
of the relative contributions of the AGN/SB components, 
because of the degeneracy introduced by the extinction of 
the AGN component: the same ratio can be obtained with 
a faint, unobscured AGN mixed with a powerful starburst, 
or with a much more luminous, even bolometrically domi- 
nant, but heavily obscured AGN. However, on average we 
observe an anticorrelation between this AGN indicator and 
the 3.3 /jm PAH starburst indicator (Fig. 3). 

The second part of our analysis consists of a quantita- 
tive estimate of the physically relevant parameters for our 
sources, i.e. the absorption of the AGN component, t, and 
the continuum fraction of emission due to the AGN, a (es- 
timated at 3.3 /im). We use the simple method described in 
R06: the total spectrum is reproduced as a combination of 
an AGN power law component, Jagn = Ci\^q N , absorbed 
by dust following the extinction law r tx A~ 1,75 (Cardelli et 
al. 1989), and a starburst component, consisting of a power 
law continuum fsB = C2^Fsb an d a Gaussian emission line 
at A = 3.3/^m with a fixed equivalent width EWsb with re- 
spect to the starburst continuum component. Following the 
results obtained for bright ULIRGs in R06, we chose the 
values Yagn = —0.5, Fsb = —0.5, and EWsb = 110 nm 
for the AGN and SB templates. 

Using this model, with simple algebra we can obtain 
the observed parameters F and EW3.3 as a function of the 
extinction r and the AGN fraction a — Ci/(Ci +C2). r and 
a can then be obtained simply inverting these equations, as 
shown in R06. The final result is: 

EWsbY - EW3.3 



TL 



' (5{EWsb - EW3.3) 

EWsb - EW3.3 

{EWsb - EW3.3) + EW3.3 * er^ 



(1) 



(2) 



where /3 is the index of the extinction power law (,3=1.75). 
The errors on a and r can be directly estimated from the 
above equations, as shown in the Appendix of R06. 

In five cases the errors on the observed parameters are 
too high to obtain significant estimates of r and a. For all 
the other sources the results are listed in Table 1. 

The application of this model to the measured pa- 
rameters listed in Table 1 is not straightforward for all 
sources: in several cases the best fit values and the 90% 
confidence intervals of one or both of EW3.3 and F are out- 
side the "allowed range" of our model, F < F m in =-0.5 and 
EW3.3 > EWmax=lW nm. In these cases our approach is 
the following: 

• If r < —0.5 at > 90% statistical significance, we assume 
t — and we estimate the relative AGN fraction using 
EW3.3: a = 1 - £W3. 3 /(110 nm); 

• If EW3.3 > 110 nm and F is not higher than the pure 
starburst value Fsb = F m in = —0.5 we assume the source is 
predominantly powered by a pure starburst and we put an 
upper limit of 1% to the AGN contribution. We note that 
in no case EW3.3 > 110 nm and F > —0.5 at a 90% confi- 
dence level. The values of a and r so obtained are shown in 
Table 1. 

The final step in our analysis is the estimate of the 
fraction of the bolometric emission due to the AGN, obol- 
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Figure 2. 3-4 (im continuum slope (in a A — fx spectrum) ver- 
sus Equivalent Width of the 3.3 fim PAH emission feature for 
the sample of z < 0.15 ULIRGs in the 1 Jy sample (Veilleux et 
al. 1999). The dashed lines mark the starburst, unobscured AGN 
and AGN zones according to the results of Risaliti et al. 2006. 
The sample is split into two different panels for clarity, depend- 
ing of the optical classification of the objects, marked as different 
symbols. 



This quantity is related to the AGN fraction in the L-band, 
a, through two new parameters, i.e. the ratio between the 
L-band and bolometric luminosities for pure AGN and pure 
SB, Ragn and Rsb- Simple calculations show that: atd = 
a /{a. + K(l — a)) where K — Ragn / Rsb ■ Following the 
results obtained for bright sources (R06) we use Ragn ~ 
0.2 and R SB = 0.002. The factor K = 100 in the above 
equation imply that the bolometric contribution of the AGN 
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Figure 3. Ratio between L-band and bolometric luminosity, Rl, 
versus equivalent width of the 3.3 /im PAH emission feature. We 
adopt 20% errors on Rl, estimated from the average IRAS and 
L-band absolute calibration uncertainties. The dashed lines mark 
the starburst and AGN zones according to R06. Symbols are the 
same as in Fig. 2. 

can be small even in those cases where the AGN component 
is dominant in the L-band. The values of ctboi estimated 
using the above equation are listed in Column 6 of Table 1. 
The uncertainties related to the factor K (log(if)~ 2±0.5, 
R06) are discussed in Section 3.3. 



3 DISCUSSION 

The application of our model to the 59 z < 0.15 ULIRGs in 
the 1 Jy sample of Veilleux et al. (1999) provided a quantita- 
tive estimate of the AGN/starburst contribution in ULIRGs, 
and in particular in those optically classified as LINERs. 

Here we (1) review and discuss our results, (2) com- 
pare them with the optical spectroscopic classification, and 
with the other multi- wavelength AGN/starburst indicators 
summarized in the Introduction, and (3) discuss the uncer- 
tainties of our approach. 

3.1 The AGN content in ULIRGs 

The most relevant result of the present work is a quantita- 
tive estimate of the contributions of the two energy sources 
to the luminosity of ULIRGs. The estimated AGN fractions 
in the L-band and in the total emission are summarized in 
the histograms in Fig. 4. The comparison between the two 
distributions shows the power of our diagnostics: the AGN, 
when present, dominates the L-band emission in most cases, 
even in objects where its contribution to the total luminos- 
ity is negligible. We note that our estimates strongly de- 
pend on the modelization of the starburst/ AGN emission. 
In particular, we assume a fixed starburst template, with no 



free parameter. We expect this to be a simplification of a 
more complex spectral behaviour, in at least two respects: 
a) the intrinsic equivalent width of the 3.3 /im emission fea- 
ture (assumed to be EW3.3 = 110 nm in our model) shows 
a significant spread among pure starbursts, and b) the ob- 
served starburst emission could be affected by some red- 
dening/absorption, as suggested by the large spread of the 
ratio between flux of the 3.3 /im feature and the bolomet- 
ric emission in pure starbursts (106). Therefore, we believe 
our analysis provides correct results in a statistical sense for 
the whole sample, but can be affected by larger uncertain- 
ties than estimated in Table 1. We will discuss this issue in 
more detail in Section 3.3. 
Overall, our analysis shows that: 

- An AGN is present in most (> 60%) ULIRGs. This frac- 
tion drops to > 30% for LINERs (however, in several objects 
with low S/N the upper limit for the AGN fraction is quite 
high, see Table 1). 

- The contribution of AGNs to the L-band emission is 
> 50%, while the fraction of the bolometric luminosity of 
ULIRGs due to AGNs is < 20%. Restricting to LINERs, 
the fraction of AGN contribution is < 30% in the L-band, 
and < 15% bolometrically. 

- In a few LINERs the presence of an heavily obscured AGN 
is bolometrically significant. On the other side there is an 
indication that most of the heavily obscured AGN are LIN- 
ERs (this is based only on ~ 5 — 8 objects, so it needs to 
be confirmed with larger samples) . This finding confirms the 
heterogeneous nature of this class of objects, which are on 
average dominated by starburst emission, but can host pow- 
erful buried AGNs. 

- In our model, QSO class luminous buried AGNs are found 
only in a small fraction of ULIRGs classified optically as 
LINERs or starbursts. 

3.2 Comparison with other AGN indicators 

An AGN/SB classification for most of our sources has been 
presented, based on observations at other wavelength, or on 
a different analysis of L-band data. 

Optical: Fig. 2 shows the classification of our sources in a 
continuum slope-PAH EW plot, for each optical class: AGN, 
starburst, and LINER. AGNs and starbursts are clearly sep- 
arated, and are located in the regions found by R06. Two ex- 
ceptions are found among optically classified starbursts: two 
sources, IRAS 20414-1651 and IRAS 11506+1331 (Panel C 
in Fig. 1) show signatures of the presence of absorbed AGNs, 
which were missed by studies at other wavelengths. LINERs 
are spread throughout the plot, confirming their composite 
nature, as quantitatively shown in Table 1 and Fig. 4. 
L-band absorption features: Aliphatic hydrocarbon ab- 
sorption features at rest-frame wavelengths 3.4-3.5 fim have 
been detected in seven sources (three AGNs and four LIN- 
ERs). All these objects have a clear detection of an AGN 
according to our model. Moreover, they are all associated to 
steep continua, with estimated continuum extinction r > 1 
(see also Sani et al. 2008). Therefore, this correlation is a 
further independent confirmation of the validity of our di- 
agnostics. 

Previous L-band diagnostics: Most of the sources in our 
sample were presented and discussed in 106. In this work, 
the possible presence of the AGN in this work is estimated 
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Figure 4. Distribution of the relative AGN contribution to the 
L-band emission (upper panel) and to the bolometric emission 
(lower panel) in our sample of ULIRGs. Blue and red (dark) his- 
tograms refer to LINERs only, while yellow (light) histograms 
refer to the whole sample. 



through the EW of the 3.3 /tm PAH and the presence of 
absorption features. The final classification is shown in Ta- 
ble 1, and is in general agreement with our results. We stress 
again the two main improvements in the analysis presented 
here: a quantitative estimate of the relative AGN/SB contri- 
butions, and the inclusion of the continuum reddening as a 
further indicator of an obscured AGN (e.g. the source C in 
Fig. 1). 

Mid-IR diagnostics: Imanishi et al. (2007) presented com- 
plete Spitzer-IRS spectroscopy of most of our sources, and an 
AGN/SB classification based on several diagnostic features 
(PAH emission at 6.2, 7.7 and 11.3 fj,m, silicate absorption at 
~9-10 ^m). The results are shown in the last column of Ta- 
ble 1, and again are in good agreement with the conclusions 
presented here, with the same limitations discussed above 
for previous L-band diagnostics. 

Nardini et al. (2008) performed an analysis of the same 
Spitzer spectra, focusing on the 5-8 /im spectral range, and 
applying an analysis similar to the one presented here. The 
results on the single sources are in excellent agreement in the 
two works: this can be verified comparing our estimates with 
their Table 1. Here we only note that the L-band analysis, 
besides providing an independent AGN/starburst deconvo- 
lution, has a higher diagnostic power, especially for ULIRGs 
hosting highly obscured AGNs, thanks to the higher contrast 
between AGN and starbursts in the L-band than at 5-8 /im. 

Alternative, powerful mid-IR indicators of the 
AGN/starburst contributions in ULIRGs are the sili- 
cate absorption features in the 9-10 fim band (Spoon et 
al. 2007), and the presence of high-ionization emission 
lines (Farrah et al. 2007). These studies, performed on 



samples with large overlaps with ours, provide results in 
good agreement with our classification. However, as for the 
other approaches described in this Section, no quantitative 
estimates have been done based on these indicators. 
X-rays: A strong X-ray emission is one of the strongest 
AGN indicators. If the obscuration of the AGN is not 
too high (column densities Nh<W 24 cm -2 ), the primary 
AGN spectrum is directly visible at energies below 10 keV, 
providing a good estimate of its total luminosity, through 
bolometric corrections obtained from AGN SEDs (e.g. 
Risaliti & Elvis 2004). In the obscuration is higher, only 
the reflection component is visible in the X-rays. In these 
cases, the spectral features of the reflection spectrum can 
identify the presence of the AGN, however a good estimate 
of its total luminosity is not possible, due to the large 
uncertainty in the reflection efficiency. Extensive studies 
of local ULIRGs (e.g. Franceschini et al. 2003, Ptak et 
al. 2003) show that the current X-ray observatories can 
provide useful spectra for this diagnostics only for the 
10-20 brighest sources, while most of the ULIRGs in our 
sample are too faint for this kind of analysis. The agreement 
between X-ray and L-band diagnostics for the few ULIRGs 
with high quality X-ray spectra is quite good (Sani et 
al. 2008, Braito et al. 2009). However, this analysis cannot 
be extended to more sources until more sensitive X-ray 
observatories are available. 



3.3 Model uncertainties 

The model presented here is based on a series of strong as- 
sumptions, which may introduce systematic errors, in addi- 
tion to the statistical errors obtained from the fitting pro- 
cedure. The main possible simplifications, already discussed 
in R06, are the assumption of fixed AGN and starburst tem- 
plates, and the absence of any absorption/reddening in the 
starburst component. Other possible relevant effects may be 
caused by different dust extinction curves from the one as- 
sumed here. 

Regarding possible obscuration of the starburst compo- 
nent, we note that, though it cannot be excluded, it is not 
strongly required by the data. Physically, this may be due 
to the origin of the starburst emission: the integrated con- 
tribution of a large number (> 10 6 ) of single sources. Each 
single star can be heavily absorbed (as observed, for exam- 
ple, in the central regions of our own Galaxy), but the single 
differences are smeared out by the large number of sources, 
producing a rather constant spectral profile of the total emis- 
sion. Regardless of the physical interpretation, this is also 
suggested by the observations: the highest quality spectra of 
pure starburst are remarkably similar, both in the L-band 
(R06) and in the 5-8 fim range (Nardini et al. 2008, Brandl 
et al. 2007). 

Regarding the fixed AGN and SB templates, this is an 
obvious over-simplification, as shown by the scatter in the 
spectral parameters (Fig. 2, and R06). We may try to esti- 
mate the effect of this scatter by changing the parameters 
of the pure AGN and pure SB templates, and reporting how 
this affects the estimate of the parameters a and r. Since 
this scatter is hard to measure in a self-consistent way, we 
prefer to adopt a different, observation-based approach. In 
order to do this, we notice that our model can easily pro- 
vide an estimate of the expected bolometric luminosity for 
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Figure 5. Upper panel: Ratio between intrinsic L-band and bolo- 
metric luminosity versus the AGN fraction a obtained from our 
analysis. The blue continuous line is that expected based on our 
model. The cross on the top left corner indicates the average 
statistical errors on the two plotted quantities. Bottom panel: 
Comparison between the estimated and measured bolometric lu- 
minosities. The 0.3 dex dispersion around the 1:1 relation is a fair 
estimate of the total uncertainty in our prediction of the AGN and 
SB bolometric contributions. 



each source, as a function of Ragn, Rsb, the model param- 
eters a and r, 
In particular, 



and the observed L-band luminosity, L®-4 ■ 



Lbol — [(xRagn + (1 — a)RsB] X L^IJ 11 



(3) 



where L^1J R is the intrinsic L-band luminosity (i.e. cor- 
rected for the extinction of the AGN component) and is re- 
lated to the observed one by: L^J H =L^_ B 4 s /[a+(l-a)e" T ]. 
Our estimate of Lbol can then be compared with the ob- 
served Lbol, as measured by IRAS, by using the standard 
calibration of Sanders & Mirabel (1996). In Fig. 5 we plot 
two relations: 

- In the upper panel the ratio Rcorr between the 
absorption-corrected L-band luminosity, L^_J R , and the 
measured bolometric luminosity is plotted versus the L-band 
AGN fraction, a. The continuous line is the expected rela- 



tion based on our model, with no free parameters. In par- 
ticular, the shape of the curve is determined by the mathe- 
matical relations between the model parameters, while the 
values at a=0 and a=l are the parameters Rsb and Ragn, 
respectively. Overall, the agreement between the expected 
curve and the data is quite good. This is a strong confirma- 
tion of the validity of our method: our simple, two-parameter 
model can successfully predict the intrinsic L-band to bolo- 
metric ratio. 

- In the lower panel we directly compare the predicted and 
measured values of Lbol- This is the most direct way to 
estimate the overall uncertainties in our analysis. These in- 
clude both the statistical errors from the spectral fits, and 
the systematic errors due to the fixed templates. The dis- 
persion from the 1:1 relation is lower than 0.3 dex. This is 
enough to provide a rough estimate of the AGN contribu- 
tion in single sources, and an accurate average estimate for 
the whole sample. 



4 CONCLUSIONS 

L-band spectral analysis of ULIRGs is a powerful way to 
disentangle the AGN and starburst contribution to the bolo- 
metric emission. This is primarily due to the diagnostic 
power of (a) the 3.3 /im PAH emission in starbursts, (b) 
the continuum slope, indicative of the reddening of the AGN 
component, (c) the high ratio (about 100) between the AGN 
and starburst emission at 3 /xm, for the same bolometric lu- 
minosity. 

We applied a simple model to the available L-band spec- 
tra of 59 ULIRGs from the 1-Jy sample, at redshift z< 0.15, 
in order to detect and quantitatively estimate the AGN con- 
tribution to the emission of these sources. We find that 
AGNs are present in 60% of ULIRGs, but in most cases 
they are not the main energy source. As a consequence, their 
overall contribution to the total ULIRG luminosity is ~20%. 
We find that the subsample of objects optically classified as 
LINERs is rather composite, with similar AGN/SB fractions 
as the whole ULIRGs sample. 

Our analysis revealed several objects optically classified 
as starburst or LINERs, and/or with weak or absent previ- 
ous evidence of the presence of an AGN in the near/mid- 
IR, which according to our model host a powerful, obscured 
AGN. We plan to observe these objects in the hard X-rays, 
in order to obtain an independent, unambiguous confirma- 
tion of these findings. 
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